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Abstract 

We study the interaction between a one dimensional magnetic nanostructure and a thin film 
superconductor. It is shown that different magnetic distributions produce characteristic magnetic 
field signatures. Moreover, the magnetic structure can induce a weak link in the superconducting 
film, or be positioned directly above a predefined nonsuperconducting weak link. We estimate the 
magnetic flux associated with such a structure, and discuss a general expression for the energy 
calculated within the London model. 
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Weak links in superconductors have generated a lot of interest over the past decades, 
both in the study of conventional and high T c superconductors. More recently, research has 
also been focused towards junctions in thin films, due to their potential in future technology 
(see e.g. Ref. QXJ and references therein). For such junctions it has been found that the 
field is a superposition of fields from Pearl vortices along the junction with a certain line 
density [I]]. An interesting special case is the generation of weak links with spatially localized 
magnetic fields. Creation of a weak link in a bulk (or thick film) superconductor using a 
magnetic domain wall was first proposed by SoninjfJ]. In that paper expressions for the 
magnetic fields where found, and it was argued that domain walls can induce movable weak 
links. In the current paper we try to extend the idea of Ref. || to thin film superconductors, 
and estimate the magnetic field and flux distribution associated with an one dimensional 
magnetic nanostructure. The weak link could be a domain wall, generated and controlled 
by a stress pattern or external field, or it could be a stationary prefabricated nanomagnetic 
stripe. A particular advantage of using magnetic domain walls as weak links is that they can 
be moved at high speeds, and may therefore have potential applications in future fluxtronics 
devices. On the other hand, in a prefabricated magnetic stripe (e.g. permalloy) the polarity 
of the magnetic vector could easily be switched by an external field 0. This could also be 
of interest for creation and annihilation of vortices. Since the vortex pinning energy often 
can be regarded as proportional to the thickness of the superconductor, it is reasonable 
to assume that pinning by nonmagnetic sources are negligible in a thin film. Therefore, 
it should in principle be possible to produce film systems in which the magnetic texture 
determines the junction properties. 

Consider a thin superconducting film located at z = with thickness much smaller than 
the penetration depth of the superconductor. The surface is covered by a one dimensional 
magnetic structure centered at x = 0, with thickness much smaller than that of the super- 
conducting film. It is assumed that the magnetic structure is sufficiently long that end effects 
are not a problem, and that it consists of surface charges separated from the superconductor 
by a very thin (negligible) oxide layer to avoid the exchange of electrons or spin (see e.g. 
Ref. 0]). We will here analyze the case where the magnetization is perpendicular to the 
plane of the superconducting film, since we expect this geometry to give a stronger coupling 
to the superconductor. In order to gain some insight into the behavior of the weak link, let 
us first compare the magnetic fields from the two following models for the magnetization 
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distribution: 

M G = M Q exp(-ax 2 )5(z)e z , (1) 

and 

M s = 



M 5(z)e z if -W < x < W 
if \x\ > W 



where M z0 and a are constants. Here we assume that a = 2W~ 2 . 

Using the method Ref . || |] , we find that the gaussian magnetization distribution results 
in the following field components (per unit length): 

rr r, , A e M f°° , 9 exp (— kl/Aa) cos(k x x) . ,,,,,, 

H?(x,z) = / k\ PK *'J h exp(-k x \z\)dk x , (2) 

q. s A e Mo f°° 7 2 exp (—kl/Aa) sin(k x x) 



H?(x,z) = '-^ ' x ' ^l^ m eM-h\z\)dk x ■ (3) 



On the other hand, the step magnetization distribution gives 

s 2A e M f°° sin{k x W)cos{k x x) 
H x (x,z) = / k x exp(-fc x z Jdfcj, , (4) 

5 2A e M Z" 00 sm(^iy)OT(ti) 

H^[x,z) = / k x — — — exp(-k x \z\)dk x . (5) 

t Jo 1 + 2A e fo c 



Figure shows iff when z = A e /200 (solid line) and z = (dash-dotted line). Note that 
when z = the magnetic field oscillates due to the steep magnetization gradient. This could 
therefore be interpreted as Gibbs oscillations, well known in Fourier analysis. At a certain 
height above the surface these oscillations are smoothed out due to the exponential decay 
factor. Figure [| shows Hf (solid line) and H x (dash-dotted line) when W = A e /40 and 
z = 0. Note that the peak of the z component is located at the origin. This is in contrast 
to Fig. U, where the maximum field is located near the edges. Moreover, the two negative 
peaks of the field are much less pronounced for a gaussian distribution. 

It is clear that if the magnetic field exceed the critical field of the superconductor, then 
a weak link is generated at which vortices may exists^. The current across the link is, in 
absence of any external currents, given by 

J s x = dJ c sin (02 - 0i + 2tt— ) , (6) 
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where = 02 — 0i + 27r$/$o is the phase difference across the junction, $ the magnetic 
flux through the junction and $o the flux quantum. The magnetic flux for a gaussian 
magnetization distribution is estimated to be 

* « „ r r HAviz = '^ M ° r kl <^tm^Mi dkx , (7) 

where Ay = y 2 — y\ is the length of the weak link. It is seen that the flux, and therefore also 
the phase, is dependent on the magnetization, the wall width and the penetration depth. 

Let us now consider a normal weak link (i.e. not magnetic), in which the phase can 
be found by using the approach of Ref. [JIJ]. If we position a very thin magnetic structure 
(as discussed above) directly over the predefined weak link, a flux will flow through this 
junction. Here we may set A e — > oo, since the junction is entirely nonsuperconducting, and 
this results in 

$ « f exp(-k 2 x /4a)sin(k x W)dk x = C/i M Ay , (8) 

^V 71 "" Jo 

where C is an unimportant constant (C = [exp(— 2)/y / 7r] exp(k x )dk x ), and we have as- 
sumed that a = 2W~ 2 . This approach also gives a reasonable description of the case where 
the strength of the magnetic field from the magnetic structure breakes down superconduc- 
tivity. However, the simple treatment given here only accounts for the direct magnetic flux 
through the junction, and neglects the vortices distributed around the weak link. In general, 
the full expression for can not be found explicitly, but can be obtained by first determining 
the vortex interspacing distance by evaluating the expression for the energy (i.e. find the 
zero energy). From the given distribution of vortices one may obtain the magnetic field H vy , 
and finally take advantage of the expression J x = dJ c sin(j) = —2H vy to obtain an expression 
for the phase 0. Such a numerical analysis is outside the scope of this Brief Communication. 
However, let us for completeness discuss some properties of the energy of the system, which 
is in general found by evaluating 

E = J^oH 2 + i/i A 2 j2 - /i M • H)dV . (9) 

Here /xo is the permeability, and the current due to supercurrents and magnetization gradi- 
ents can be written as 

J = J s + J m = J s + V x M (10) 
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Note that the integration over the energy density is taken over the whole space, although the 
current can only flow in the volume of the thin film superconductor. I.e., we do not adopt the 
usual approach of dividing the space into superconducting and nonsuperconducting regions. 
Thus, we believe that the approach shown below can be applied to a more general class of 
systems, as long as there is no exchange of electrons and spin between the magnetic and 
superconducting structures. In order to obtain a more useful expression for the energy, we 
first transform the part associated with kinetic energy of the superconducting electrons 

E kineUc = J \ 2 AdV = J \ 2 J S - [V x (H - M)] dV . (11) 

Here we have used that J s = V x (H — M). To further transform this integral, we note 
that a surface integral over the kernel J s x (H — M) vanishes when the surface is located 
far from the system. This means that we can write 

E ki netic = \^ J \ 2 J 2 s dV = J y X2 (H- M) V X J s dV . (12) 
But now it should be remembered that the London equation gives 



vxJ s = -^ + ^y, (13) 



where V is the vortex source function, which represents all the vortices in the system (In 
the case of Pearl vortices it is simply a sum of delta functions). We obtain 

E ki netic = J \ 2 AdV = J [-H 2 + M-H + (H-M)-V]dV . (14) 

In total, the energy becomes 

E = -hq [ [-M ■ H + (H - M)-V\dV . (15) 
2 Jv 

Equation IT51 can be applied to systems which are sufficiently local that the surface correc- 



tions can be neglected. The usefulness of Eq. [15] relies on the fact that for a thin film 
superconductor it reduces to a two dimensional integral, since the vortex source function 
and the magnetization distribution both are assumed to be located at z=0. Then the Fourier 
analysis of Refs. [|], || can be applied to obtain simple integrals for the energy. In their 
interesting paper Erdin and coworkers demonstrated a different method for calculating the 
energy within the London approximation 0]. The technique presented here gives an alterna- 
tive and perhaps more intuitive route to evaluate the energy. It is instructive to divide the 



energy terms in three different parts. The self-energy of the magnetic structure is 

E n = ~fiojM-H m dV , (16) 
the energy associated with the vortex is 

E v = ^fMiJ H v -VdV , (17) 

and the interaction energy 

E vm = ~fi J^[M-H V - H m V + M'V] dV . (18) 

It is seen that the E vm consists of the interaction between the magnetization and the vortex 
field, the magnetically generated field and the vortex source function and also the inter- 
action between the magnetization and the vortex source function. Note that the two first 
contributions have opposite sign, and their sum is equal to —fiQ j v M • H v dV only in the 
case when they are equal. Although this is a reasonable approximation in some particular 
cases, the above analysis shows that it is not correct in general |5|, 0, ||. Moreover, one should 
take into account the finite distribution of the vortex source function. 

In conclusion, we have discussed some properties of thin film superconductors in the 
close vicinity of one dimensional magnetic structures. It was found that the magnetic field 
depends strongly on the magnetic distribution function, and simple estimates for the flux 
through a weak link were derived. Finally, we derived an expression for the energy associated 
with a more general class of systems. 
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FIG. 1: The z component of the magnetic field generated by a step-like magnetization distribution 
when z = (dash-dotted line) and z = A e /200 (solid line). Here W = A/40, and the curves have 
been normalized with respect to the maximum peak of the dash-dotted line. 
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FIG. 2: The x (dash-dotted line) and z (solid line) component of the magnetic field generated 
by a gaussian magnetization distribution. Here z=0, W = A/40 and a = 2W~ 2 . The curves are 
normalized with respect to the maximum peak of the z component. 



